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Spiral model, jamming percolation and glass-jamming transitions 
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Abstract. The Spiral Model (SM) corresponds to a new class of kinetically constrained models introduced 
in joint works with D.S. Fisher [8,9]. They provide the first example of finite dimensional models with 
an ideal glass-jamming transition. This is due to an underlying jamming percolation transition which has 
unconventional features: it is discontinuous (i.e. the percolating cluster is compact at the transition) and 
the typical size of the clusters diverges faster than any power law, leading to a Vogel-Fulcher-like divergence 
of the relaxation time. Here we present a detailed physical analysis of SM, see [5] for rigorous proofs. 
We also show that our arguments for SM does not need any modification contrary to recent claims of Jeng 
and Schwarz [10]. 
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1 Introduction that a particle does not (or cannot) move if surrounded by 

too many others. This can be also interpreted in terms of 

Theoretical progress in understanding- the glass and jam- , . r -i-i ^- r^i ^^^ t^i^t^jt u j- u • 

^ ° t> t> J dynamic lacnitation [4]. All KCM snare two basics prop- 

ming transition, and more eenerally glassy dynamics, is i- i / • \ / n- \ i -r 

° o J o J J erties: particles (or spins) can move (or nip) only ii a 

hampered by the shortage of finite dimensional models , . , • ; ; i i r • i • i i 

certain constraint on the number oi occupied ncignbors 

that display the basic phenomenological ingredients and . .^i,^ .i -i^ t ,^ i -i 

IS verined. Unce the constraint is verified tire dynamical 

that are simple enough to be fully analyzed. Kinetically , i ^ i^- n ii i- ^ -i ^■ 

" J J J rules are such that the resulting hsoltzmann distribution 

Constrained Models (KCM) [1] are an exception. They i ^ j r ^ a 

^ ' I- ' ^ -'is trivial. I.e. uncorrelated from site to site. As a conse- 

have been introduced few decades ago [2.3] as models for ii i ^ -i- -r • i i • i 

<- ■ ' quence the glass transition, it any, is purely dynamical 

glass-forming liquids. They are based on the assumption 
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in these models. Furthermore, another advantage is that real glass or jamming transition should display (if they ex- 
the study of the dynamical transition can be reduced to a ist) . They have been also rigorously proved in [5] modulo 
highly correlated percolation problem (of a new kind). In the standard conjecture on the existence of two different 
fact, in these models a particle (or a spin) can be blocked correlation lengths for DP [11]. 

if it has too many occupied (up) neighbors which can be In the following we will sketch in a physical and hope- 
blocked by their neighbors and so on and so forth. Using fully transparent way the arguments which lead to the 
the fact that the Boltzmann equilibrium distribution is above results providing the tools needed to analyze this 
trivial one can prove [6,7] that the only dynamical transi- transition and explaining the underlying mechanism: it is 
tion that can take place in these systems corresponds to a the consequence of two perpendicular directed percolation 
jamming percolation where an infinite cluster of mutually processes which together can form a compact network of 
blocked particles (spins) appears. In [8, 9] we introduced frozen directed paths at criticality. In the final section we 
a new class of KCM which displays such a transition on shall discuss the generality of our approach and the uni- 
a finite dimensional lattice, we will thus refer to these versality of the jamming percolation transition of SM. 
models as Jamming Percolation (JP) models. Here we will 
review the easiest example of a JP model, namely the two- 
dimensional spin model which has been introduced in [5,9] 
and dubbed Spiral Model (SM). For SM the existence of a 
jamming transition has been rigorously proved [5] and the 
exact value of Pc has been identified: pc coincides with the 
critical threshold of directed site percolation (DP) in two 
dimensions, p^^ ~ 0.705. Contrary to recent claims [10] 
our proof for SM does not need any modification (we will 
pinpoint in section [4.11 the incorrect assumption of [10]). 
This jamming transition has remarkable properties: the 
density of the frozen cluster, <?(p), is discontinuous at the 
transition but the cross-over length over which the sys- 
tem is still ergodic (or liquid) diverges. Furthermore the 
time scale for relaxation (and also the cross-over length) 
diverges faster than any power law. These properties are 
quite unusual but are exactly what is often assumed the 



2 The spiral model and its related percolation 
problem 



Consider a square lattice and, for each site x, define among 
its first and second neighbours the couples of its North- 
East (NE), South- West (SW), North- West (NW) and South- 
East(SE) neighbours as in Fig.l, namely NE= (a:; + ei,a;-|- 
61 + 62), SW= (x— 62, 61 — 62), NW= (a;— 61, X— 61 -I- 62) 
and SE= (2; + 6i, a; + 6i — 62), where 61 and 62 are the coor- 
dinate unit vectors. The spiral model is a stochastic spin 
lattice model where a spin can flip if and only if the fol- 
lowing constraint is verified: both its NE and/or both its 
SW neighbours are down and both its SE and/or both its 
NW neighbours are down too (see Fig.l). If the constraint 
is verified then the spin flip rate is p from down to up and 
1 — p from up to down. As a consequence, the invariant 



Giulio Biroli, Cristina Toninelli: Spiral model, jamming percolation and glass-jamming transitions 



NE a) NE 



b) 3 Critical density 



NW 



o 

• o 



o o o 

NW 

o] o 

SE 



o 



SE 3.1 Occurrence of blocked clusters for p> 
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In order to establish pc < I we will identify a set of 

Fig. 1. Site X and its NE, NW, SE and SW neighbours. Filled 

blocked clusters and show that they exist with probability 

(empty) dots stand for up (down) spins. The constraint is (is 

one (with respect to the Bernoulli measure) for p > p^^ , 

not) verified at x in case (a) (in case (b)). 

therefore Pc < p^^ < 1. 

Let us start by recalling the definition and a few basic 

results on DP (see e.g. [11]). Take a square lattice with 
probability measure reached for p < pc is the Bernoulli randomly (independent) occupied sites and put two ar- 
product measure, i.e. independent from site to site. It is rows going out from each site x towards its neighbours in 
such that a spin is up with probabihty p and down with the positive coordinate directions, x+ei and x+e-z- On this 
probabihty 1 - p. As explained in the introduction the dy- directed lattice a continuous percolation transition occurs 
namical transition takes place when an infinite cluster of at a non trivial critical density pf-^ ~ 0.705 (a percolating 
mutually blocked spins appears with probability one with cluster is now one which spans the lattice following the di- 
respect to the Bernoulh product measure. An easy way rection of the arrows). This transition is second order, as 
to unveil the existence of this cluster is to run the ccUu- for site percolation, but belongs to a different universality 
lar automaton defined by the following update rule: down class. In particular, due the anisotropy of the lattice, the 
(up) spins are mapped to empty (occupied) sites; empty typical sizes of the incipient percolating cluster in the par- 
sites remain empty, occupied sites get emptied only if the allel (ei -I- 62) and transverse (ei - 62) directions diverge 
kinetic constrain is verified. The remaining cluster coin- with different exponents, ^|| ~ (p^^ -p)""^" and ^|| ~ 
cides with all spins up that are mutually blocked under with ~ 1.74 and z ~ 1.58. 

the stochastic dynamics. Note that the kinetic rules can Back to the Spiral Model, let us consider the directed 
be also rephrased by saying that at least one among the lattice that is obtained from the square lattice putting 
four sets NEUSE, SEUSW, SWUNW and NWUNE two arrows from each site towards its NE neighbours, as 
should be completely empty. From this perspective (and in Fig. 2 a). This lattice is equivalent to the one of DP, 
identifying an occupied site with high density regions in simply tilted and squeezed. Therefore, for p > p^^ , there 
a liquid) SM encodes in a very simplified way the cage exists a cluster of occupied sites which spans the lattice 
effect that emerges in liquids and granular media close to following the direction of the arrows (cluster inside the 
the glass and jamming transition. continuous line in Fig. 2 a)). We denote by NE-SW clusters 
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the occupied sets which follow the arrows of such lattice 
and NW-SW clusters those that follow instead the arrows 
drawn starting from each site towards its NW neighbours. 
Consider now a site in the interior of a spanning NE-SW 
cluster, e.g. site x in the Fig. 2 a): by definition there is 
at least one occupied site in both its NE and SW neigh- 
bouring couples, therefore x is occupied and blocked with 
respect to the updating rule of SM. Thus, the presence of 
the DP cluster implies a blocked cluster and pc < p^^ 



9 O • 




o o o 



a) 




b) 



Fig. 2. a) The directed lattice obtained drawing arrows from 
each site towards its NE neighbours. Particles inside the con- 
tinuous line belong to a NE-SW spanning cluster, thus they 
are blocked, b) A NE-SW non spanning cluster crossing two 
NW-SE clusters. We depict both types of T-junctions: the bot- 



foUows. Note that these results would remain true also for tom (top) one occurs with (without) a site in common for the 

a different updating rule with the milder requirement that NE-SW and NW-SE cluster. In both cases the last point of 

only at least one among the two couples of NE and SW NE-SW before the crossing is blocked by the presence of the 

sites is completely empty (and no requirement on the NW- NW-SE cluster, thus all sites inside the dashed line cannot be 

SE direction). However, as we shall see, the coexistence of erased before erasing at least one site of the NE-SW cluster, 
the constraint in the NE-SW and NW-SE directions is cru- 



cial to find a discontinuous transition for SM, otherwise 
we would have a standard DP-like continuous transition. 



a T-junction with NW-SE paths, as shown in Fig. 2 b) (see 
section |4T1 for a detailed definition of T-junction). By us- 
ing such T-junctions it is also possible to construct frozen 
clusters which do not contain a percolating DP cluster nei- 
ther in the NE-SW nor in the NW-SE direction: aU NE- 
SW (NW-SE) clusters are finite and arc blocked at both 
Before showing that below p^^ blocked clusters do not oc- ends by T-junctions with finite NW-SE (NE-SW) ones (see 



3.2 Absence of blocked clusters for p <p\ 



DP 



cur, a few remarks are in order. If instead of SM we were Fig. 3 b)). As we will show in section [4?2l these T-junctions 

considering the milder rules described at the end of previ- are crucial to make the behavior of the transition for SM 

ous section, the result would follow immediately since the very different from DP, although they share the same crit- 

presence of a blocked cluster would imply the existence of ical density. This also means that the fact that spanning 

a DP one. On the other hand for SM rules, since blocking DP clusters do not occur for p < pf^ is not sufficient to 

can occur along cither the NE-SW or the NW-SE direction conclude that also blocked clusters do not occur. What 

(or both), a directed path implies a blocked cluster but the strategy could one use? Recalling bootstrap percolation 

converse is not true. This is because a NE-SW non span- results [12], a possible idea is to search for proper unsta- 



ning cluster can be blocked if both its ends are blocked by ble voids from which we can itcratively empty the whole 
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lattice. Of course, since we already know that blocked clus- 
ters occur when p > p^^ , something should prevent this 
unstable voids to expand at high density. 





Fig. 3. a) Qi (continuous line). Si (dashed line) and the four 
sites (*) which should be emptied to expand Qe. We draw 
the necessary condition for x to be frozen: it should belong 
to a NE-SW cluster spanning Si and supported by NW-SE 
cluster from the exterior. b)The frozen structure described in 
the text: continuous lines stand for occupied NE-SW or SW- 
NE clusters. Each of these clusters is blocked since it ends in a 
T-junction with a cluster along the transverse direction. The 
dotted rectangle adjacent to cluster AB (EF) are the regions in 
which this cluster can be displaced and yet a frozen backbone 
is preserved: the T-junctions in C and D (G and H) will be 
displaced but not disrupted. 



Consider the region inside the continuous line in 
Fig. 3 a), namely a "square" of linear side \/2l tilted of 45 
degrees with respect to the coordinate axis and with each 
of the four vertexes composed by two sites. If Qe is empty 
and the four sites external and adjacent to each vertex 
denoted by * in Fig. 3 a) are also empty, then it is possible 
to enlarge the empty region Qe to Qe+i- Indeed, as can 
be directly checked, all the sites external to the top right 



side can be subsequently emptied starting from the top 
one and going downwards. For the sites external to the 
other three sides of Qe we can proceed analogously, some 
care is only required in deciding whether to start from top 
sites and go downwards or bottom ones and go upwards. 
Therefore we can expand Qe of one step provided all the 
four * sites are empty or can be emptied after some iter- 
ations of the cellular automaton. Let us focus on one of 
these * sites, e.g. the left one, in Fig. 3 a). As it can be 
proved by an iterative procedure (see [5]), in order for xl 
not to be emptyable there should exist a NE-SW cluster 
which spans the square Se of size £ containing the top left 
part of Qe (region inside the dashed line in Fig. 3 a)). This 
is due to the fact that any directed path in the NW-SE 
direction can be unblocked starting from the empty part 
of Se below the diagonal in the ei + 62 direction. There- 
fore the only way in which x^ can be blocked is that it 
belongs to a NE-SW cluster that is either supported by 
NW-SE clusters running outside Se or that is infinite. In 
any case this NE-SW cluster has to be at least of length 
^. As a consequence, for large £, {£ » ) the cost for 
a one step expansion of is proportional to the proba- 
bility of not finding such a DP path for any of the four * 
sites, 1 — 4exp(— c£/^||). Thanks to the positive correlation 
among events at different £'s, the probability that the emp- 
tying procedure can be continued up to infinity is bounded 
from below by the product of these single step probabil- 
ities which goes to a strictly positive value for p < p^^ 
since ^|| < 00. Note that, as we already knew from the re- 
sults of section 3.1, this is not true for p > p^^: the pres- 
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ence of long DP paths prevents the expansion of voids. As blocked, as long as the NW-SE paths are blocked. If the 

a conclusion, the probability of emptying the whole lat- T-junction occurs with a site in common (site inside the 

ticc starting from an empty square Qi centered around a square of Fig. 3 b) this is a trivial consequence of the fact 

given point in the lattice and with I » ^|| is finite (al- that this point belongs to the NW-SE path. In the other 

though very small). Since there is an infinite number of case it can be easily checked that the last point belonging 

points in the lattice, there will be at least one (actually a to the NE-SW path (site inside the circle of Fig. 2 b) is 

finite fraction) of sites from which the whole lattice can blocked thanks to the one above it, which belongs to the 
fl 

be emptiecu for p < p^^. This, together with the result NW-SE path. All other possible crossings are related to 
of section 3.1, yields pc = p^^ . these two cases by symmetry. In [10] it is stated that in 

order for a NE-SW path to stabilize a NW-SE path (or the 
converse) they shouldn't only cross but also have a point 

4 Critical behavior 

in common and since this may not happen our proof for 
4.1 T-junctions the SM needs a modification. As explained above this con- 



clusion is incorrect and our proof for the SM model does 
not need any modification (see also [5] for further details). 



4.2 SM: Discontinuity of the transition 



One of the most important characteristic of the SM model, 
already alluded to in the previous sections, is that a di- 
rected path in the NE-SW direction can be supported by 
another path running in the NW-SE direction (and vicev- 

ersa) via a T-junction. Let us discuss this point in detail In the previous sections we have shown that the perco- 
since recently it has been incorrectly claimed that this lation transition due to the occurrence of a frozen back- 
is not true for the SM model. There are only two pos- bone for the Spiral Model occurs at pf^. We will now 
sible types of crossing of a NE-SW path with a NW-SE explain why the density of the frozen cluster is discon- 
one: either they have one point in common or not. In the tinuous, 'Pip^^) > (the frozen structures are compact 
latter case they should cross as in the upper crossing of rather than fractal at criticality) . 

Fig.2 b). In both cases we call the crossing a T-junction By translation invariance <P{p^^) is equal the proba- 

and the key observation is that if a NE-SW path ends in bihty that a given point, e.g. the origin, is blocked (i.e. 

two T-junctions with NW-SE paths (or viceversa), it does it belongs to an infinite blocked structure). In order to 

not need to continue beyond the crossings in order to be show that <P(pf^) > we will then construct a set of 



^ Mathematically, one would say that the ergodic theorem configurations, B, for which the origin is blocked and such 
implies that in the thermodynamic limit with probability one that Pts{p^^) > 0. Since <P{p^^) > Pig our result implies 
the final configuration is completely empty. <P{p^^) > 0. In order to define B, consider a configu- 
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ration in which the origin belongs to a NE-SW path of 
finite length £o/2: this occurs with some finite probability 
qo > 0. Now focus on the infinite sequence of pairs of rect- 
angles of increasing size £iX £i/12 with £i = £q, £i ~ 2£i^2 

, . , . • ^ \ A n -11 Fie. 4. a) The sequence of intersecting rectangles, b) Dotted 

and intersecting as m iig.4 a). A connguration belongs 

non straight line stand for NE-SW (NW-SE) clusters spanning 

to B if each of these rectangles with long side along the 

the rectangles c) Frozen structure containing the origin 

NE-SW (NW-SE) diagonal contains a NE-SW (NW-SE) 




percolating path (dotted lines in Fig. 4 b)). If this is the , . ; • , • r i i i • i n i . 

^ o // analyzing the typical size, Lc, below which frozen clusters 

case then the infinite backbone of particles containing the „ . , , , , . 

occur on nnite lattices. 

origin (cluster inside the continuous line in Fig. 4 c)) sur- 

We first obtain a lower bound on Lc constructing ex- 

vives thanks to the T-junctions among paths in intcrsect- 

plicitly blocked structures that exist with finite probability 



ing rectangles. Therefore <P{p) > qo Yli^i ^ P{£i)'^, where 

length s intersecting as in Fig. 3 b). This type of structure 



as long as i < L[^. Consider NE-SW and NW-SE paths of 

P{£i) is the probability that a rectangle of size £{ x 1/I2£.i 



with short side in the transverse direction is spanned by ^ 

can be emptied completely only starting from its border 

a DP cluster. Recall that there is a parallel and a trans- 

since each finite directed path terminates on both ends 

verse length for DP with different exponents, i.e. a cluster 

into T-junctions with a path in the transverse direction. 

of parallel length £ has typically transverse length £^ [11]. 

Therefore it is frozen if we continue the construction up 

Let us divide the £i x I/I2£i rectangle into £l~^ slices of 

to the border of the lattice and we take periodic bound- 
size £i X I/12£|. For each slice the probability of having 

ary conditions. Thus the probability that there exists a 

a DP cluster along the parallel direction at p^^ is order 

frozen cluster, I — R(L,p), is bounded from below by the 

unity. Thus, the probability of not having a DP cluster 

probability that each of the 0{L/s) dotted rectangles in 

in each of the slice is I — P{£i) ~ 0[exp(— c^|~^)]. From 

Fig. 3 b) contains at least one path connecting its short 

this result and the above inequality we get ^(pf^) > 0. 

sides. This leads to R{L,p) < (L/s) exp(— cs ^) pro- 

Therefore the infinite cluster of jamming percolation is 

vided s < (for s > C|| the probability of having a DP 

"compact" with dimension = 2 at the transition. 

cluster in a rectangle starts to decrease and cannot be 
bounded anymore by I — 0[exp(— c£,- ~^)]). Thus taking 
4.3 Crossover length s oc C|| , we get log >ki\p- pf ^ | . 

In order to establish an upper bound on Lc, we deter- 
Let us now focus on the divergence of the incipient blocked mine the size L"^ above which unstable voids, that can 
cluster approaching the transition. This can be studied be expanded until emptying the whole lattice, occur typi- 
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cally. The results in Section [3H imply that the probability 5 DisCUSSion 

of expanding an empty nucleus to infinity is dominated 

Let us first discuss the dynamical behavior of the SM 

by the probability of expanding it up to £ = ^|| . Indeed, 

model. The results of the previous sections have impor- 

above this size the probability of an event which prevents 

tant consequence on the dynamics of SM. First, incipient 

expansion is exponentially suppressed. Therefore, consid- 

blocked clusters can be unblocked only from the boundary. 

cring the 0(L/^||)^ possible positions for the region that 

As a consequence the relaxation timescale is expected to 

it is guaranteed to be emptyable up to size ^|| , we can 

scale at least as (but likely larger than) their typical size: 

bound the probability that a lattice of linear size L is 

(X exp(A:/|p — p^^\^)- Indeed this can be proved rigor- 

emptyable as R{L,p) > L^S, where 6 is the probability 

ously [7]. Furthermore, since the fraction of blocked sites 

that a small empty nucleus can be expanded until size 

is finite at the transition, two point dynamical correlation 

^11 . In the emptying procedure described in Section [3. II we 

functions, e.g. spin-spin correlations, will show a plateau 

evaluated the cost for expanding of one step the empty re- 
like super-cooled liquids approaching the glass transition. 

gion Qi. Analogously, the cost of expanding directly from 

The plateau, also called Edwards- Anderson parameter in 

Qe to Q2i can be bounded from below by , with 

the context of spin-glasses, corresponds to the frozen spin 

C a positive constant independent from £. This can be 

fluctuations. These two dynamical characteristics are re- 
done by dividing the region contained in Q2i and not in 

markable since they lead to a dynamics qualitatively sim- 

Q( into £^^'' strips with parallel and transverse length 

ilar to the one of glass-forming liquids. It would be very 

of order i and requiring that none of them contains 

interesting to perform more detailed investigations and 

a DP path which percolates in the transverse direction 

comparisons, in particular by numerical simulations. 

and using for this event the scaling hypothesis of directed 

The extension and universality of the jamming per- 

percolation when p p^^ . Thus for the expansion up 

eolation transition of SM remain fundamental questions 

to size we get 6 > [iS^" c^'*'"^' ~ exp(-C'e|j""), 

to be investigated. As it has been discussed in [9] (see 

with C" > . This, together with above inequality, yields 

also [10, 13]), it is possible to identify a class of rules which 

logLc<fc«b-pf^r''. 

give rise to a jamming transition and belong to the same 
universality class of SM: as p Pc the divergence of the 
incipient frozen cluster follows the same scaling and the 
As a consequence upper and lower bound leads to the transition remain discontinuous. For all these models the 
same scaling at leading order and imply that the cross- jamming transition is a consequence of the existence of 
over length diverges with an essential singularity, i.e. faster (at least) two transverse directed percolation (DP)-like 
than any power law for p p^^ . processes which can form a network of finite DP-clusters 
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blocked by T-junctions with clusters in the transverse di- 
rection. A model that belongs to this class is for example 
the Knight model defined in [8]. Note that in general, at 
variance with SM, it will not be possible to determine an- 
alytically the exact value of pc (this was possible for SM 
thanks to the fact that in each of the two transverse di- 
rections the blocked clusters can be exactly mapped into 
those of conventional 2 dimensional DP). Neither it is pos- 
sible to generalize the rigorous proofs obtained for the SM. 
However, it is nevertheless possible to obtain numerically 
a reliable estimate of pc and a confirmation that the tran- 
sition has the same properties of SM. This is done analyz- 
ing finite size effects with proper choices of the geometry 
and boundary conditions which allow to focus separately 
on each of the two transverse directions. In this way one 
can verify that on long length scales the two indepen- 
dent directional processes are in the DP universality class 
and obtain a good numerical estimate of pc- Using suit- 
able boundary conditions and geometries is particularly 
important for jamming percolation since, as for bootstrap 
percolation, convergence to the asymptotic results can be 
extremely slow. For an extended discussion on this we refer 
to [9] , where the value of the critical density for the Knight 
model has been derived. The result is p^'^'^aht ^ 0.635 
and differs from our original conjecture p^™sht _ pDP jgj 
which was due to the overlooking of some blocked struc- 
tures [13]. 



iming percolation and glass-jamming transitions 9 

References 

1. F.Ritort, P.SoUich Adv. in Phys. 52 (2003), 219. 

2. G.H. Frederickson and H.C. Andersen, Phys. Rev. Lett. 53, 
1244 (1984). 

3. W.Kob, H.C.Andersen, Phys. Rev. E 48 (1993) 4364. 

4. J. P. Garrahan and D. Chandler, PNAS 100, 9710 (2003). 

5. C. Toninelli and G. Biroli, Spiral Model: a new cellu- 
lar automaton with a discontinuous glass transition arXiv: 
0709.0378 

6. G. Toninelli and G. Biroli. J. Stat. Phys. 117, 27 (2004). 

7. N. Cancrini, F. Martinelli, G. Roberto, C. Toninelli Kineti- 
cally constrained spin models, to appear in Probab. Th. and 
Rel. Fields, preprint [math.PR/0610l(j6j 

8. G. Toninelli, G.Biroli, D.S.Fisher, Phys. Rev. Lett 96, 
035702 (2006) 

9. G. Toninelli, G.Biroh, D.S.Fisher, Phys. Rev. Lett. 98, 
129602 (2007). 

10. M. Jeng, J. M. Schwarz , On the study of jamming perco- 
lation, larXiv:0708.0582l 

11. H.Hinrichsen, Adv. in Phys, 49 (2000), 815; R. Durrett, 
Ann. Prob. 12, 999-1040. 

12. R.H.Schonmann, Ann. of Probab. ,20 (1992), 174-193; 

13. M. Jeng, J. M. Schwarz,Pi7S. Rev. Lett. 98, 129601 
(2007). 



We thank D.S.Fisher for very important collaborations 
on this subject. 



